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Inflationary models involving more than one scalar field naturally produce isocurvature pertur- 
bations. However, while these are fairly well studied, less is known about their evolution through 
the reheating epoch, when the inflationary flelds decay into the standard constituents of the present 
universe. In this paper, by modelling reheating perturbatively, we calculate the power spectrum 
of the non-adiabatic pressure perturbation in three different inflationary models. We show that 
the isocurvature can grow large initially, but decays faster than the pressure perturbations. When 
reheating ends, the isocurvature is negligible for the double quadratic and double quartic inflation- 
ary models. For the product exponential potential, which features large isocurvature at the end of 
inflation, the isocurvature decays during reheating and is around five orders of magnitudes smaller 
than the pressure perturbation at the end of reheating. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The inflationary paradigm is incredibly successful in 
solving the original problems of big bang cosmology [I] , 
to such an extent that it is now firmly a part of standard 
cosmology. However, despite this success, we have no sin- 
gle, specific theory of inflation which fits the data better 
than all others. There has been a plethora of inflation- 
ary models considered in the literature to date, ranging 
from the original and simplest single field, slow-roll po- 
tential [SHI] , through more general models involving mul- 
tiple scalar fields (see, e.g., [S] and references within), to 
more complicated scenarios involving scalar fields with 
non-standard kinetic terms, such as k-infiation [H]. It 
is anticipated that data obtained from current and fu- 
ture experiments, such as observations of the cosmic mi- 
crowave background (CMB) from the Planck satellite[7], 
will enable us to constrain the inflationary paradigm. 

One signature which can be used to distinguish be- 
tween different inflationary models is the isocurvature, 
generated by the relative entropy perturbation between 
scalar fields. This naturally occurs in any general system 
containing more than one field. Constraints on primor- 
dial isocurvature come from observations of the CMB. 
However, to date, the data only constrain the ratio of 
isocurvature to adiabatic perturbations to be of the or- 
der of 0.1 [5]. An interesting possibility for determining 
a favoured inflationary model from its primordial isocur- 
vature signature is through the vorticity induced by non- 
adiabatic pressure at second order in perturbation the- 
ory [5141 1|. The possibility of constraining isocurvature 
in this manner is discussed in more detail in Ref. [X2\. 

In addition to the unknowns around the specific model 
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of inflation, we do not know how the scalar fields driving 
inflation decay in the standard model particles that exist 
in the universe today. This process is dubbed 'reheat- 
ing' and there is currently no agreed upon mechanism by 
which the universe reheats. 

For models consisting of multiple scalar fields (exclud- 
ing those where the isocurvature has decayed by the end 
of infiation since these can then be treated as single field 
systems) this process of reheating is of utmost impor- 
tance. This is because isocurvature perturbations nat- 
urally 'spoil' the standard picture of the early universe 
since the curvature perturbation is no longer conserved 
on super-horizon scales. This makes it interesting, and 
important, to address how the isocurvature remaining at 
the end of the inflationary phase evolves through a period 
of reheating. 

There are different methods of modelling the reheating 
phase, either by focussing on the particle interactions and 
particle production during reheating (see, e.g., Ref. ^3] 
and references therein), or by considering, in a more clas- 
sical manner, how the inflationary fields decay into the 
matter and radiation components of the standard cosmo- 
logical model. It is the latter technique which we adopt 
in this article, using full cosmological perturbation the- 
ory p^HTS] . Of course, if we consider the scalar fields 
decaying into a single fluid, as has been done in some 
other work, the system will not be able to support any 
evolution of isocurvature. 

Therefore, in this article we consider the decay of the 
inflationary scalar fields into more than one species of 
fiuid and investigate how reheating affects isocurvature 
perturbations generated during multi-field inflation. Us- 
ing results from Ref. ^12,, we model the reheating pro- 
cess by a system containing two scalar fields along with 
a matter and radiation fiuids with decay constants con- 
verting the scalar fields into fluids. By modelling the sys- 
tem like this, we explore the scenario in a more general 
manner than in previous work, and consider the non- 
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adiabatic pressure perturbations between the radiation 
and the matter fluids after reheating is complete. In do- 
ing so, we can get a handle on whether it is common for 
isocurvature to survive through reheating. 

This paper is organised as follows: in the next sec- 
tion, we review the details of perturbative reheating, and 
present the set of equations to be solved. In Section [lIT| we 
review the numerical procedure. We present the results 
for the three different models in Section ITVl and conclude 
in Section IVl 



II. THE MODEL 

In this paper we consider linear perturbations to a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time and work in the flat gauge, in which the line element 
takes the form 

ds^ = -{l + 2(f>)dt'^ + 2aB^^dx'dt + a'^{t)6^jdx'dx^ . (2.1) 

Here (f) is the lapse function, and B is the scalar shear; 
we consider flat spatial slices in agreement with current 
observations [S]. Our index convention is as follows: 
{/i, ■ ■ ■} run over 4 space-time indices, {i, j, . . .} denote 
the 3 spatial indices, {a, /?,...} are labels for the 'fluids' 
(including matter and radiation and the kinetic and po- 
tential fluids attributed to the scalar fields - see later) 
and {A, B, . . .} label the scalar fields. 

In this work we model reheating perturbatively, fol- 
lowing the elementary theory of reheating developed in 
Refs. [inilin]- This amounts to including an additional 
friction term, -^Tip, in the equation of motion of the scalar 
field during reheating, i.e.^ 



^+{3H+-T)cp + U,, 



0. 



(2.2) 



However, this approximation is only valid when the field 
is rapidly oscillating about its minimum, and the addi- 
tional friction term should not be present during infla- 
tion. In our analysis, we ensure that the approximation 
is valid by only allowing the decay constants for a par- 
ticular field to be non-zero once that field has passed 
through its potential minimum and is taking part in re- 
heating (see Sec. HI). We describe the system in more 
detail in the following. 

In order to model the system through the reheating 
phase, we consider a system of multiple fluids, where the 
total energy-momentum tensor, 7),,^, containing the fields 
and the fluids, is covariantly conserved, V^T'^i/ = 0. The 



^ Note that our decay constants differ from those in, e.g., Refs. |21l 
I22| by a factor of 1/2. This difference arises due to the introduc- 
tion of the decay constants in the fluid conservation equations as 
opposed to the field evolution equations. 



individual fluids, on the other hand, have a modified con- 
servation equation, 

(2.3) 



where the energy-momentum transfer functions sum to 
zero: qI""* = . The energy-momentum transfer vec- 
tor has the components [El [23l [24] 

Qio.)o = -QM + ^)-SQa, (2.4) 

= {fa + QaV)^, . (2.5) 

We now need to specify a functional form for the var- 
ious transfer components. Let us consider a system con- 
taining N scalar fields, with a Lagrangian 

ji: = IY.^a + u{^b), (2.6) 

^ A 

along with pressureless matter and radiation fiuids, m, 7, 
respectively. We adopt the approach of Ref. [53] and split 
the scalar fields into kinetic fiuids, with energy density 
and pressure pA = Pa = i and a single potential 
fiuid with pu = U = —Pij . We can then define the energy 
transfer functions, following, e.g., Ref. [53], which can be 
expressed explicitly in terms of field variables as 

QA = -\{Tt + T^y\-U^A^A. (2.7) 

Qu^Y. ^^^^^ ' (2.8) 

A 

^ A 

Q-^ = \Y.^^^I- (2-10) 



Here, the V are decay parameters (i.e. constants, apart 
from r^, which is proportional to the potential of the 
scalar fields). 

The energy transfer functions can then be perturbed 
to give 

U,A' 



5Q, 



ri + r^ + ^){^ASvA-^l^ 



(PA 

'PA^U^A^Sipj , 



W = 



X! X! ^-Aj^^j + "^U.AiSlfA - 



5Q^ 



A 

YV^iipAS^PA 



Va<4 



fA^i 



(2.11) 



(2.12) 
(2.13) 

(2.14) 



Finally, we take the /„ to be constants (and hence the 
fa.i vanish)^, in agreement with, e.g., Ref. [53]. 



In fact, we can only take the /'s transferring momentum from 
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A. Background 

We now consider the evolution of the various fluid com- 
ponents. In the background, energy conservation gives 
one equation for each fluid, a, 



Pa + 3H{pa +Pa) =Qa 



(2.15) 



So, the radiation and matter fluid conservation equations 
are 



P7 



^ A 



(2.16) 
(2.17) 



where we use the equations of state for the radiation and 
matter: — \p-^-, Pm — 0. The fields each have a 
Klein-Gordon equation coming from energy-momentum 
conservation of the kinetic fluids 



^a + Ua^O, (2.18) 



and the evolution equation for the potential fluid is sat- 
isfied identically. 



where the evolution of Sq^n and Sq-y is governed by their 
respective Eq. (2.201 and, for the fields, we have Va = 



—Slpa/^a giving the expression 



V 



PmVm + ^P-fV^ - J2a VA^Vj 



Pm + 3P7 + "Pa 
Then, the fluid energy conservation equations are 



(2.23) 



6p^ + AHSp^ 



4V^ 
3 a 



-VjP^ 



1 



'Pa<P 



(2.24) 



Spm + SHdpm H WmPn 

a 



J2'^m('PASiPA - ^^a4' 
A ^ ■ 

(2.25) 



and momentum conservation gives 



5q-y + 3HSq^ + -p^(i 



Sqm + 3HSq 

n 



Sp-y V 



(2.26) 



V 



(2.27) 



B. Linear perturbations 

Considering now the evolution of linear perturbations, 
the energy-momentum conservation equations yield the 
governing equations for each (perfect) fluid, a, [24], 

Spa + 3H{Spa + SPa) H Wq (pc -|- Pa) 

a 

= Qa^ + 5Qa, (2.19) 
5qa + iHSqa + {pa + Pc.)'t> + - Qc^V) = , (2.20) 

where the momentum perturbation, 5qa — a{va + 
B){pa + Pa) and Va is the three- velocity of the a fluid. 
V is defined as 



V 



For our system this becomes 

V = 



Sqm + Sqy + J2a VaVa 



Pn 



3^7 + La Va 



(2.21) 



(2.22) 



the fields into the fluids to be zero - due to the split of the flelds 
into kinetic and potential parts, the /'s controlling the transfer of 
momentum from the kinetic to the potential fluid representations 
of the flelds are non-zero. However, since we do not need to use 
the flelds' Euler equations in this work, we omit the detail here. 



In order to obtain the set of equations in a closed form 
we choose the base variables for the fluids to be the mo- 
mentum and energy density perturbations, 5qa and 5pa- 
We then need to invoke the perturbed Einstein field equa- 
tions which give the two required equations |24j . 



V2 

H — B 
a 



3H'U + H—B = -AnGSp , 
H(j) = -4ttG{p + P)V, 



(2.28) 
(2.29) 



where the total energy density perturbation is Sp — 
^a^Pa- Considering the matter conservation equation, 
Eq. ( |2.25 1, we can use the fact that Vm = V,n/a — B, to 
write 



V2 V2 
6pn^ + iHSpra H ^<5(?m ^Pn 

= ^^n^{;PA5vA - ^'Pa<I> 



Then, from Eq. (2.281, we have 

4nG 



-B 



H 



and so 



'iHSp 



V2 



6p - 3H(j) , 



(2.30) 



(2.31) 



^Qm + ^^PmSp + 3H(l)p„ 

J2ri[^ASvA~l^A^): (2.32) 
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where the total energy density is [T7| 
5p = (5pm + (5p^ + ^ (^ipA5vA-'P\(l> + U^A5(pAj ■ (2.33) 

A 

We obtain a similar equation for the radiation fluid 



Sp^ + AHSp. 



IBttG 

•2 



PjSp + AH(f)pj 
Y,V^(^^A5^A-\^A<t)- (2.34) 



Finally, to close the system we can write the lapse func- 
tion in terms of fluid and field variables, using Eq. ( 2.29 1, 
as 



47rG/ 



(2.35) 



Evolution equations for the scalar field perturbations 
can be obtained using the energy conservation equation 
with the energy density of each field given by 5pA = 
tfA5^A'~'f\4> ■ We obtain a modified Klein-Gordon equa- 
tion for the LfA field: 



5^/ 



3i/ + i(rf 



2U, 



U^AsSifiB 

3 

: V2 



"PAfl 



-vaVa =0. 
(2.36) 



Then, using the field equations the Klein-Gordon equa- 
tion for each field tpA becomes 



SifA + 

+ ^ U,ab5(Pb 
AttG 



5Vi 



B 

(pA<j>+ ^^(PaSp = 0. 



2Ua + ^MK 



H 



where, as above, the lapse function, 
by using Eq. 



-5ipA 

- r^) + iHipA 

(2.37) 

I, can be replaced 



The closed set of equations to be solved is therefore 



Eqs. (2.34), (2.32), (2.26), (2.27) for the fluids, and 



Eqs. (2.37), for the fields, along with the expression for 



from Eq. (2.35). 



C. Non-adiabatic perturbations 

For the remainder of the paper, we consider the spe- 
cial case of a system containing two scalar fields ip and 
X, though we note that the expressions can be readily 
extended to the multiple field case. 



In a general fluid system, one can define the non- 
adiabatic pressure perturbation by expanding the pres- 
sure, as 



SPwad = SP - cjp , 



(2.38) 



where the adiabatic sound speed is defined as = P j p. 
The initial non-adiabatic pressure perturbation in this 
system containing two scalar fields is given in terms of 
the field variables by [12] 



SttG 



nad 



3^2 v-.'^^ --:x 
2([/^^<5^ + U.,,8x) 

m (^2 + ^2) 



(2.39) 



<y9(5(^ + x^x + U,^^^ + U,x^X 



This can then be compared to the pressure perturbation 
which, for a system of two fields, is 

8P = + x^x - {^^ + X^)4> ~ U,^6^ ~ U^^Sx . 

(2.40) 

When the fields have decayed into the fluids, after re- 
heating, the non-adiabatic pressure perturbation is 



SPnad =^3*^/^7 



1 - 



1 pj6p„ 



Pm + Pj 

where we have used the fact that 



3 Pm + Pj 



E 



Pa 2 
P " 



(2.41) 



(2.42) 



This is the non-adiabatic pressure perturbation which we 
use to calculate the entropy perturbation remaining af- 
ter reheating for different potentials in the next section. 
Throughout the reheating phase, when the fields are de- 
caying into the fiuids, the non-adiabatic pressure pertur- 
bation will have components from both Eqs. (2.39) and 



III. NUMERICAL METHOD 

The results in Section |IV| were derived using the 
Pyflation numerical package [57]. The implementa- 
tion of the inflationary calculation is described in detail 
in Refs. [121128]. 

The Pyflation package has been updated to include 
the evolution equations of radiation and matter fiuids 
at background and first order in perturbations. These 
fiuids initially are set to zero and there is no transfer 
from the infiationary fields until each field has reached 
its minimum value for the first time. 

Reheating is assumed to end when the energy density 
of the scalar fields is reduced to below some specified 
fraction of the total energy density of the system. This 
value can be set as a parameter (piimit in the following) 
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and for the results below is either 10~^ or 10~^°. The 
scalar fields are assumed to play no further role in the 
evolution after this point and are set to zero from then 
onwards. This allows the code to avoid wasting time fol- 
lowing extremely small oscillations of the negligible scalar 
fields. The perturbations of the scalar fields are also set 
to zero to maintain the consistency of the perturbative 
scheme. Computations of the non-adiabatic pressure per- 
turbation and other derived quantities now include the 
radiation and matter fluids as described above. 



IV. RESULTS 

In this section we present our results. We focus on 
the three inflationary potentials we studied in Ref. [H], 
namely the double quadratic, double quartic and prod- 
uct exponential models. The decay constants are crucial 
to our result, since they control how quickly the fields 
decay into the radiation and matter. We will focus on 
constraints for these parameters for each specific model 
below. However, one constraint common to all models is 
that, in order for matter to be subdominant until matter- 
radiation equality, we require [211 129] 



< 



T 



decay 



< 10" 



decay 



(4.1) 



In the following we therefore ensure that this bound is 
satisfied by taking = 10~®r^ for each scalar field, A. 



A. Double quadratic inflation 

The first model we consider is the two-field generali- 
sation of quadratic inflation, the double quadratic infla- 
tionary model, with the potential 



1 



2 2 



(4.2) 



We use the parameter values from Ref. [H], choosing 
= Irriip and setting = 1.395 x 10~^Mpl. In Fig.[l] 
we plot the power spectra of the pressure perturbation, 
5P (red solid line) and the non-adiabatic pressure pertur- 
bation, (5Pnad (green dashed line). In all graphs in this 
section results are plotted against number of efoldings af- 
ter the end of inflation (i.e. inflation ends at the origin on 
the a;-axis). We can see that, in this model, as inflation 
ends the spectrum of the non-adiabatic pressure pertur- 
bation is many order of magnitude smaller than that of 
the pressure perturbation, as presented in Ref. ,12J . 

In Fig. [2] we zoom in to the last few efolds of the previ- 
ous plot, and consider only the evolution through reheat- 
ing. We can see that during the initial phases of reheat- 
ing, the isocurvature'^ perturbation sharply increases, due 
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1: A comparison between the power spectra of 5P 
(red solid line) and (5Pnad (green dashed line) for the double 
quadratic potential at the Wmap pivot scale. 




FIG. 2: A comparison between the power spectra of 5P 
(red solid line) and JPnad (green dashed line) for the double 
quadratic potential during the reheating phase only. 



to the generation of radiation and matter fluids and the 
corresponding fluctuation of the adiabatic sound speed 
of the total fluid (as seen in Fig. [S]). 

However, the non-adiabatic pressure perturbations de- 
cay faster than the pressure perturbation during the pe- 
riod of reheating. When the fields are turned off and 
reheating is over, approximately 2.5 efolds after the end 
of inflation, the non-adiabatic pressure perturbation is 
around twenty orders of magnitude smaller than the pres- 
sure perturbation. Thus, we conclude that the reheating 



Note that, in this paper, we wiU use the terms 'isocurvature' and 



'non-adiabatic pressure' interchangeably. 
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FIG. 3: The adiabatic sound speed of the system for the 
double quadratic inflation model. 



FIG. 5: A comparison of the (5Pnad power spectrum for two 
values of pumit, the value of the field energy density below 
which the field is turned off. 




FIG. 4: The (5Pnad power spectrum for different values of the 
decay constant. 



phase does not change the prediction that the isocur- 
vature is negligible in the double quadratic inflationary 
model. 

In order to model reheating for this potential we have 
had to specify the decay constants. The results in Figs, 
[l] and [2] are obtained by choosing the decay constants 
' ~' ' ' ' - - - ^ - j-^^ ^ Allowing 
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we can 



and to be equal and set to 7 x 10 ^Mp 
these to be different from one another and setting 
1 X lO-'^MpL while keeping = 7 x 10~'^AfpL 
see from Fig. [4] that this simply means that reheating 
takes longer, while the spectrum obtained at the end is 
the same magnitude. 



Furthermore, as described in Section III in the numer- 



ical procedure we drop a field from the dynamics when 
its energy density becomes smaller than a certain frac- 
tion of the total energy density. In Fig. [5] we plot the 



non-adiabatic pressure power spectrum for two values of 
this ratio, piimit = 10^^,10^^", and we can see that al- 
lowing the fields to take part in the simulation for longer, 
i.e. choosing a smaller value for piimit, results in a non- 
adiabatic pressure perturbation spectrum with a smaller 
amplitude. 

Finally, we can confirm that the Universe reheats prop- 
erly in this model by considering the complete back- 
ground evolution. Plotted in Fig. [6] are the fractional 
energy densities of the constituent components of the 
cosmic fluid compared to the total energy density. We 
can see that as the model reheats, the proportion of en- 
ergy density in the fields drops to zero and is converted 
into radiation and a small amount of matter. At this 
point in time, the adiabatic sound speed of the system is 
Cg — 1/3, as shown in Fig. [Sj as the system enters a ra- 
diation dominated era, followed by a matter dominated 
era towards the present day. We do not consider the evo- 
lution of the perturbations for the entire history of the 
Universe since our approximations are not valid, and we 
do not include dark energy perturbations. However, ob- 
taining the correct background evolution into the matter 
dominated era is a good confirmation that the model is 
sufficiently accurate. 



B. Double quartic inflation 

The next model we consider is a special case of hybrid 
inflation, the double quartic model whose potential takes 
the form 



(4.3) 
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FIG. 6: The evolution of the energy density of the fields (red FIG. 8: A comparison between the power spectra of SP (red 
solid line), radiation (green dashed line) and matter (blue solid line) and i5Pnad (green dashed line) zoomed in to the 
dotted line) as a fraction of the total energy density. reheating phase of the double quartic inflationary potential. 
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FIG. 7: A comparison between the power spectra of SP (red 
solid line) and SP^ad (green dashed line) for the double quartic 
infiationary potential at the Wmap pivot scale. 



where v,^ and ipc are constants. In order to compare 
with our previous results, we use the same parameter 
values as in Ref. ,12j . 

Plotted in Fig. [7] is the comparison between the power 
spectra of the pressure and non-adiabatic pressure per- 
turbations as a function of time. The difference from 
the double quadratic potential during the inflationary 
epoch is apparent, and we can see the sharp increase of 
the isocurvature towards the end of inflation. When fo- 
cussing on the reheating phase, as presented in Fig.|8j we 
can see that, while the non-adiabatic pressure spectrum is 
within a few orders of magnitude of that of the pressure 
perturbation at the end of inflation, as reheating pro- 
gresses, the isocurvature decays more rapidly. When this 



model has fully reheated, we can see that the isocurva- 
ture spectrum is again negligible, roughly twenty orders 
of magnitude smaller than the spectrum of the pressure 
perturbation. 

The decay parameters used for this model follow 
Ref. [21] closely, and are F!^ = = 3 x 10-"MpL and 
= = 3 X We have checked that, as in 

the previous case, varying these parameter values does 
not drastically change the results. 

In a similar way to in the previous model, we can show 
that the reheating of the double quartic model also pro- 
duces the correct background evolution of the universe, 
i.e. the fields reheat resulting in a radiation dominated 
epoch followed by matter domination. 

C. Product exponential 

Finally, we study the product exponential potential, 
a model popular in the literature due to its tendency 
to source large non-gaussianities and which has a large 
isocurvature perturbation present at the end of inflation. 
The potential is 

V = Vo^^e-^^" , (4.4) 

and we use parameter values as in Ref. P]. In Fig. [9] 
we show the results for this model, plotting the power 
spectrum of the pressure perturbation compared with the 
power spectrum of the non-adiabatic pressure perturba- 
tion. In this figure we can see the characteristic growth of 
the non-adiabatic pressure perturbation power spectrum, 
peaking around the end of inflation where it is then larger 
than that of the pressure perturbation as inflation ends. 
This was shown in Ref. [H]. Choosing the decay con- 
stants used in Ref. [21], we investigate the evolution of 
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FIG. 9: The power spectra of 5P (red solid line) and 5Pnad 
(green dashed line) for the product exponential potential at 
the Wmap pivot wavenumber. 
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FIG. 10: A comparison between the power spectra during 
reheating of 5P (red solid line) and SP^ad (green dashed line) 
for the product exponential potential. 



this isocurvature through reheating. 

In Fig. [To] we compare the pressure and non-adiabatic 
pressure perturbation for this model, focussing on the pe- 
riod of reheating. We see that after inflation ends both 
power spectra decay. The isocurvature spectrum is larger 
than that of the pressure perturbation immediately after 
inflation, but decays faster. When reheating ends the 
spectra stop oscillating and enter a period of monotonic 
decay. The non-adiabatic pressure spectrum is around 
five orders of magnitude smaller than the spectrum of 
pressure perturbations at this time, and continues to de- 
cay faster than the pressure perturbation spectrum. 



V. DISCUSSION 

In this paper we have studied the period of reheat- 
ing, where the scalar fields driving inflation decay into 
the standard matter and radiation fluids in the universe 
today, for three different two-field infiationary models. 
Each model contained a canonical kinetic term along with 
a specific form for the potential. We modelled the reheat- 
ing phase perturbatively, introducing constants which al- 
lowed for the decay of the two scalar fields into both 
matter and radiation, and solved the system fully nu- 
merically. Doing so enabled us to describe the evolution 
of the non-adiabatic pressure perturbations beyond the 
inflationary phase, and through reheating, thereby ex- 
tending work previously undertaken in Ref . |12j . 

We found that the isocurvature perturbations are sub- 
dominant in all cases. In the double quadratic model, 
the isocurvature is negligible after inflation, and after the 
reheating period remains negligible. Therefore, the pre- 
diction of this model is robust, namely that there is no 
isocurvature for a double quadratic potential. The double 
quartic potential is a model in which the non-adiabatic 
pressure spectrum is within a few orders of magnitude of 
the pressure perturbation spectrum when inflation ends. 
During reheating, the isocurvature rapidly decays, and 
when the fields have decayed into fluids, the isocurvature 
spectrum is negligible, around twenty orders of magni- 
tude smaller. The product exponential potential is the 
only model we studied in which the non-adiabatic pres- 
sure power spectrum is larger than that of the pressure 
perturbation for a short time when inflation ends. As 
the universe reheats, the isocurvature again decays faster 
than the pressure. However, when the reheating phase is 
over, the isocurvature spectrum is within four orders of 
magnitude of the pressure perturbation spectrum. 

The product exponential model is the only potential 
studied here in which the isocurvature fraction after re- 
heating is non-negligible. This differs from the prediction 
directly from inflation, however shows that isocurvature 
perturbations remain important in this case. 

An important requirement of the model of reheating 
is that it reproduces the correct expansion history of the 
universe. Thus, we ensured that this was true for the 
models studied in this paper (see, e.g.. Fig. [6]). The fact 
that the isocurvature is subdominant likely arises from 
the constraint on the decay parameters ensuring that the 
radiation/matter fraction is correct in order for big bang 
nucleosynthesis to take place. While our model gives 
us the correct energy density fractions for the history 
of the universe, we do not trust its predictions too far 
into the radiation dominated era, since in order to con- 
sider the period around recombination, we need to use a 
CMB Boltzmann code. However, this is a good check of 
the model. Of course, non-adiabatic pressure perturba- 
tions arising from the inflationary /reheating mechanism 
are not the only way in which isocurvature can present 
itself in the universe. At a later time, around matter- 
radiation equality, this relative entropy between matter 
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and radiation sources an isocurvature perturbation, as 
studied in Ref. [30] . 

In summary, in this paper we have addressed the prob- 
lem of how to obtain reaUstic predictions for inflation- 
ary models which exhibit a large amount of isocurva- 
ture when inflation ends, presenting results for the non- 
adiabatic pressure perturbation that exists when the uni- 
verse has reheated. Such isocurvature could have im- 
portant observational consequences, such as through the 
generation of vorticity at second order in perturbation 
theory (see discussion in Ref. [T^] and Refs. [SHIT]). This 
will be explored in a future publication. 
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